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Abstract 

Aperture synthesis radio telescopes generate images of celestial bod- 
ies from data obtained from several radio antennas. Placement of these 
antennas has always been a source of interesting problems. Often, sev- 
eral potentially contradictory objectives like good image quality and low 
infra-structural cost have to be satisfied simultaneously. 

In this paper, we propose a general Minimum Variance Method that 
focuses on obtaining good images in the presence of limiting situations. 
We show its versatility and goodness in three different situations: (a) Plac- 
ing the antennas on the ground to get a target Gaussian UV distribution 
(b) Staggering the construction of a telescope in the event of staggered 
budgets and (c) Whenever available, using the mobility of antennas to 
obtain a high degree of fault tolerance. 

1 Introduction 

The technique of interferometric aperture synthesis has brought about a revolu- 
tion in the field of radio astronomy in the past four decades [T2]. The technique 
works by synthesizing images from signals obtained from antennas spread over 
a large distance. For a quick introduction to basic interferometric radio astron- 
omy, please see the Appendix A. 

Most radio telescopes today use the concept of interferometric aperture syn- 
thesis. The Very Large Array (VLA) in the USA with 27 antennas and the 
Giant Meter wave Radio Telescope (GMRT) in India with 30 antennas are good 
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examples of such telescopes. A multi-purpose new generation radio telescope 
called the Square Kilometer Array (SKA) has been proposed (see Section [2] and 
the website http://www.skatelescope.org), and is currently in the design phase. 

An important problem in the field of interferometric radio-astronomy is 
to find the antenna placement on the ground that generates a required UV- 
distributionQ This problem is a computationally difficult one [10]. In addition 
to antenna placement on purely scientific merit, other dimensions are added 
when we consider logistical and financial issues. For example, given a particular 
placement, what is the optimum cable layout? What trade-off can be struck 
between scientific merit and wire length and other infrastructure costs? 

In this paper we look at three important problems that arise chiefly from 
the standpoint of astronomic merit: 

• The first is a general problem. It turns out that a Gaussian distribution 
of UV points in the radial direction and a uniform distribution in the az- 
imuthal direction is desirable for Gaussian beams [2] . It has been observed 
empirically that a uniform random placement of antennas on the XY plane 
yields a radially tapering UV distribution [TS] . A simple calculation shows 
that this is as expected. However, it may not be sufficiently close, par- 
ticularly when some antenna positions are not feasible (say, because of 
geopolitical constraints). How does one rectify this? 

• The second is motivated primarily by financial considerations. Suppose 
the funding for the project is spread over a period of time as is likely to be 
in the case of the SKA. It is also reasonable to assume that the construction 
of the entire array will take time, because of logistical reasons. In such a 
scenario, one would still like to perform experiments and get good quality 
images in the meantime. Can we come up with a placement schedule that 
enables such a graceful upgrading? 

This question has been raised in the past in this context [TO] and oth- 
ers 1 14] . A more recent version of this question concerns expanding existing 
arrays [8 . Their solution involves evaluating every antenna configuration. 
This is expensive in the scenario where the number of antennas is large, 
and there is limited time for rearrangement of antennas between successive 
experiments. 

• The third concerns mobile antennas. While there are a limited number M 
of antennas available, there are several more pads that can host an antenna 
each. Ideally, an experiment would require antennas on at least some of 
those pads. Consider the case when only some N of those antennas are 
available (due to maintenance reasons or involvement with other experi- 
ments, the rest are unavailable). On which TV pads out of the required 
pads should we place the antennas? 

We propose a Minimum Variance Method (MVM) that tackles the above 

1 For more on the terminology and concepts, please see Appendix lAl 
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problems. It is interesting to note that the same framework is useful for these 
seemingly different problems. 

Informally, this technique involves choosing N out of M (where M > N) 
possible locations for placing the antennas. To achieve this, we start by placing 
an antenna on each of the M locations. Then, we iteratively remove M — N 
antennas, one at a time, such that we stay "close" to the target UV distribution. 
Our solution takes 0((M — N) ■ M 3 ) time, which is an improvement over brute 
force solutions. A brute force solution would involve comparing distributions 
that result from all ( A ^f_ N ) choices for antenna removal. 

We conduct our experiments on a random placement scheme, on the place- 
ment data of 120 antennas/stations of the proposed Australian SKA (see next 
section) , and on the placement data of the 24 antennas of the Sub-millimeter 
Array at Mauna Kea, Hawaii (see Section [4| . 

Some relevant literature is discussed in the next section. In Section [3] we 
present the Minimum Variance Method. Section [4] shows how the technique 
applies to the three problems stated above, along with the experimental results. 
Section [5] concludes the paper. 

2 Previous Work 

This work has been done keeping in mind the proposed Square Kilometre Array 
(SKA) project. The SKA is an ambitious multi-purpose new generation radio 
telescope with a total collecting area of approximately one square kilometer, 
designed to work over a wide range of radio frequencies - 70 MHz to 25 GHz. The 
telescope is expected to play a major role in answering key questions in modern 
astrophysics and cosmology. This high resolution array will be 50 times more 
sensitive, and will be able to cover the sky 10000 times faster than any imaging 
radio telescope array previously built. It is hoped that besides exploratory 
astronomy, the SKA will provide insights into many interesting questions about 
the birth and evolution of galaxies, origins of magnetism, possibility of existence 
of life on some other planet, verification of general relativity etc. For more details 
see g]. 

A lot of work has been done in the past on the problem of antenna placement. 
They range from constructing the array ab-initio to incrementing an existing 
one. We cite a few important example papers here; the reader is encouraged to 
read the citations therein. 

In 1989, Treloar compared various array configurations in terms of the 
amount of sampling of different regions in the UV-plane and the absence of 
holes therein, for fixed declinations [T5]. Among the configurations studied was 
the spiral, with higher concentration towards the centre with distribution of an- 
tennas varying with the radius. Cohanim et al posed the array design problem as 
a multi-objective optimization problem of maximizing image performance and 
minimizing cable- length using genetic algorithms and simulated annealing [10 . 
Motivated by the work of Takeuchi et al [Hj, they mention the problem of 
phased deployment of the antennas. 



3 



Lonsdale and Cappallo make the case for a large number of antennas, of 
the order of a thousand [9 . Among other factors including high fidelity and 
good UV-coverage, they argue that a large number of stations (a collection 
of antenna imaging a specific region of the sky) relieves the dependence on 
the earth's rotation for a comprehensive sampling of the UV-plane. They also 
present log-spiral and hybrid (log-spiral for inner regions and pseudo-random 
star for the outer regions) antenna layouts, which have the advantages of good 
UV-coverage towards the centre, a wide range of baseline lengths and shorter 
communication cables in general. Problems of insufficient coverage and long 
cable lengths appear in the outer regions. Several variants of a technique that 
seeks to maximize the distance between UV points have also been explored in 
the past (pQ and 0). 

Karastergiou et al [8] adapted the approach of Boone [1] and Cornwell [6] to 
the case when there are more potential sites (called pads) than antennas, and 
one has to choose a configuration for an experiment by shifting the antennas 
among the pads. The work stemmed from an observation of Boone [2J that if 
the density of the UV-plane is sparse, it is advisable to spread the Gaussian 
along a radial direction a bit so that all regions of the UV-plane have at least a 
sample point. The technique they used was inspired by the physical phenomenon 
of charges spreading on a closed 2-D surface in order to minimize the energy. 
Their idea is to affect all the discrete shifts and see which configuration results 
in the least energy, thus satisfying the experimental requirements. 

Su et al [T3] proposed the uniform weight approach. The problem that they 
tackle is of choosing sites for placement of antennas given several more candidate 
sites. They aim for a uniform and complete UV-distribution in the absence of a- 
priori information on the objectives of the experiments. Their method involves 
assigning a weight to a UV point that is equal to its distance from the nearest 
UV point. An antenna is thus assigned the sum of all distances of the relevant 
UV points. The antenna that has the least weight is dropped. 

Finally, we mention the tomographic projection method of de Villiers [7] in 
which he compares the projection of an imperfect and an ideal UV distributions 
onto one dimension. The discrepancy yields correction terms that are mapped 
to new antenna positions. Doing this in several directions allows one to get close 
to the ideal configuration. 

3 The Minimum Variance Method 

The basic idea of this paper is as follows. First, we divide the UV-plane into p 
regions. Division is necessary because we wish a distributed removal of points 
so that some portion of the UV plane is not denuded of UV points. The general 
way to capture the notion of distribution is to grid the plane and then talk of 
"points per grid unit" 0H3]- How this division is done will be described towards 
the end of this section. We note, however, that the MVM can accommodate any 
scheme for division of the UV plane in principle. The desired UV distribution 
dictates the number of UV points each region should hold. 
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Our algorithm starts with M > N antennas and removes M — N antennas, 
one at a time, in an iterative fashion. Suppose that at the start of the r th 
iteration, there are p^ regions and antennas. Then, JlfM — 1 UV points 
will be dropped as a result of the removal of one antenna during the r th iteration. 

Depending on the requirements of the UV distribution, we would like the 
removal of R\ UV points from the i th region in the r th iteration. However, in 
general, the number of UV points dropped from the region i due to removal of 
antenna j would be different, say Wi(j). 

We capture the discrepancy between the actual and ideal UV points dropped 
in every region by removing the antenna j with the following term: 

Var{j)=Y,^i{j)-R\ ] f- 

We now drop the antenna that has the least Var{.) value. In short, the logic 
for deletion of the N antennas is this: 



1. while there are more than N antennas remaining do 

2. for each remaining antenna j 

3. Calculate Var(j). 

4. end for 

5. Remove the antenna that has the least Var(.). 

6. end while 



This antenna removal routine can be tweaked for different problems. Let us 
now estimate the worst case time complexity of this routine. There are at most 
M — N antennas to be removed. Every remaining antenna is a candidate for 
removal, of which there are at most M. For every such candidate we go through 
the (at most M 2 ) UV points, identify those associated with the candidate, and 
determine which region they belong to, in order to calculate Var(j). Thus, we 
need at most 0((M - N)M(M 2 + p)) = 0((M - N) ■ M 3 ) steps. 

But for a rotation for changing coordinates, a UV point is defined as u = 
Xl ~ X2 and v = Yl ^ Y2 for every pair {X\,Yi) and (X 2 ,Y 2 ) of XY coordinates 
(see Appendix A). Therefore it is sufficient to work with the distribution that 
the XY coordinate differences may follow. In all our experiments, without loss 
of generality, we take A = 1 . 

There are at least two natural choices for the shape of the regions into which 
the UV plane can be divided. One is into concentric circles due to radial sym- 
metry. The other is into coaxial ellipses because of the facts that the UV tracks 
are elliptical in general, and that in several earth rotation aperture synthesis 
telescopes, a prolonged coverage is desirable. In a scheme that removes antennas 
progressively from a large set, there are two options regarding the number of 
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regions. The first is to adaptively reduce the number of regions in accordance 
with the reduced antennas. In doing so, the configuration in an iteration will 
mimic that in the previous iteration. The second option is to stick to the same 
number of regions throughout, so that all intermediate configurations will mimic 
the final. We demonstrate the use of different UV plane division schemes, both 
in terms and shape and number of regions, in our experiments. 

4 Applications 

4.1 Problem 1: UV to XY 

The UV distribution that is desired in most cases is a Gaussian along the radial 
and uniform along the azimuthal direction. As the next lemma indicates, a 
random distribution of antennas on the ground, that is, the XY-plane, gives a 
tapering distribution in the UV-plane. 

Lemma: Consider a S x S square grid having (S + l) 2 lattice points in 
the XY plane. Suppose that we choose each lattice point independently with 
probability < p < 1 and look at the UV distribution generated by the chosen 
points. Then the expected number of UV points at the position (x,y), where 
-S < x,y < S, is P \S+ 1 - \x\)(S + l - \y\). 

Proof: For the given S x S square grid, the number of pairs of lattice 
points (Xi, Y\) and (X 2 , Y 2 ) such that X 2 — X x = x and Y% — Y\ = y is (S + 
1 — \x\)(S + 1 — \y\) . Since each such contributing pair is expected to be chosen 
with probability p 2 , linearity of expectation gives the desired result. □ 

The above lemma indicates that a random XY distribution induces more UV 
points corresponding to short XY distance pairs. This tapers off for longer dis- 
tances. With this as a starting point, we apply the Minimum Variance Method 
to arrive at the target Gaussian. 

Given: (i) A geography, that is, dimensions of the land on which to place 
antennas, (ii) The number of antennas that are to be installed and (iii) The 
amplitude a/ of the desired Gaussian UV plane. 

Goal: To find a placement of N antennas that yields a Gaussian UV distribution 
centered at the origin having an amplitude a,f and standard deviation a = B/A, 
where B is the largest inter UV point distance)^] 

Approach: We start with an overkill of antennas. Instead of the required N, we 
start with M « 2N antennas and place them randomly on the XY-plane. Our 
random choice results in a distribution of points on the UV-plane that deviates 
from the desired Gaussian. We now try to minimize the deviation from the 
target Gaussian by iteratively removing an antenna that causes the minimum 
deviation from the desired Gaussian. We stop when we have removed M — N 
antennas. 

2 We work with a = B/A, since the area under the gaussian distribution curve before the 
4er coordinate closely approximates the total area. 
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Figure 1: The steeper curve is generated by an initial random distribution of 
the M antennas while the pink curve is the desired Gaussian. 



Figure [T] illustrates conceptually the line of approach. The two meshes in- 
dicate the starting and the final distributions. The effort will be towards min- 
imizing the discrepancy between the distributions at every iteration when we 
remove an antenna. 

Details: As stated earlier, our random choice results in a distribution of points 
on the UV-plane that deviates from the desired Gaussian. 

We then draw a fixed number p = B/d of concentric circles with the radii 
increasing in arithmetic progression of common difference d. We will show how 
to fix d shortly. Given these annuli, it is possible to plot a density histogram of 
the number qj 1 " 11 ' of UV-points that lie in the annuli bounded by the circles i and 
i + 1 (for 1 < i < p; gj* mt ' being the number of points in the central circle) versus 
distance. Thus, the total number of UV points initially is Q' init ' = 53i=i ?!'"''' — 
M ■ (M — l)/2 and area under this histogram is A^ intt ^ = d ■ ^i=i Q 1 ]""'' — 
d - M ■ (M -l)/2. 

Let us now fix d. The area under half the Gaussian centered at the origin 
with amplitude a and having a standard deviation of a is given by aay^ (see 
Appendix B). Thus, given a/, a and the final number N of antennas, we have 

rw _ d-N-(N-l) T , r , _ a / g %/¥ _ a / B ^f 

u / u V 2 — 2 ' J-neieiuiCj u, — jy.(jy-i) — 2JV-(JV-1) ' 

As has been mentioned previously, the central idea behind the present ap- 
proach is to crop the initial curve to fit the Gaussian. Let the number of antennas 
remaining at the end of the r th iteration (1 < r < M — N) be MM = M — r. 
Then, M^l = M. Note that removal of an antenna during the r th iteration re- 
sults in the removal of M^^ 1 ^ — 1 UV-points, which in turn results in a reduction 
of the area under the curve from A^'^ to = d ■ M M • (AfM - l)/2. 

We try to minimize the deviation from Gaussian by iteratively removing the 
antenna that causes the maximum deviation from the desired Gaussian. We use 
AW - d ■ MM • (MM - l)/ 2 = aoJ\ , to get a = HM-r). ( M-r-i)/2 



2 ' 

~ V 2 

With a determined (and a already being fixed), the expected Gaussian is 
uniquely defined. The vertical coordinate in the Gaussian distribution at x = i 

_d-(i-l/2) 2 [r— 11 

is given by the formula ae 2^ . If q\ is the number of UV points 
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remaining in region i at the end of iteration r — 1, then the number of UV 
points we would like to remove from i th region in r th iteration is 



R 



fr _H d-(M-r)-(M-r-l)/2 
l\ e 



For each region i, let Wi(j) be the number of UV-points that would be deleted 
from the region i if the antenna j is removed. Then, as per the Minimum 
Variance Method, we choose that antenna j for removal which has got the least 
value of 

/ '\ , \r-i] d-(M-r)-(M-r-l)/2 *-v-im' 
Var(j) = Y J (w l (j)-W 1 — -e )) 2 . 

For the experiment, we set a goal of generating a Gaussian distribution on 
the UV plane using N — 28 antennas with a,f = 70. We start with an initial set 
of M = 48 antennas that are generated randomly as shown in Figure 2(a)| The 
geography that we use, and the randomly generated antenna positions yields 
B sa 11180 wavelengths. Then, d ss 6470 wavelengths. 

Figure 2(c) shows the UV-distribution generated by these antennas. The 



Density Histogram of the distribution is shown in Figure 2(e) Notice that the 



histogram deviates from the desired Gaussian at several places. 

Figure |2(f)| is the cropped histogram arrived at finally. Observe that this 
histogram is closer to Gaussian than the initial one. Figures [2(d) | and |2 (b) | show 
the final UV and XY distributions respectively. 



(a) The initial set of antennas 



(b) The final set of antennas 



4.2 Problem 2: Staggered Construction 

Statement: A telescope of M antennas has been proposed with all the sites 
identified. However, the telescope construction has to be staggered and it has 
to be built in phases since the funding becomes available in small chunks over 
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(c) The initial UV-distribution 
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(d) The final UV-distribution 



258 




(e) The initial UV-density histogram (f) The final UV-density histogram 

Figure 2: UV to XY 



a time period. Suppose that we are given a budget of N antennas for Phase 1. 
Which N of the total M antennas should we construct in Phase 1? 

Ideally one would like to start making quality observations after completion 
of Phase 1 itself which means that the UV-distribution that is generated by 
Phase 1 antennas should be a good approximation of the final UV-distribution 
that we would get after placing all N antennas. 

If M is even, then we choose p^ = M — 1 ellipses with M/2 UV-points per 
ellipse. Else, we choose p^ = M ellipses with (M — l)/2 UV-points per ellipse. 
One ellipse is chosen for every M of all the elliptical UV-tracks sorted by their 

9 2 

u + si ^ 2 g . We aim for the removal of the same number of UV points from each 
region during every iteration. Thus, = M p[!-f 1 , and 

w 

p mW - 1 

= £(«*(?) - — w- )2 - 

i=l " 
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4.2.1 Experimental Results 



We ran the experiment on the scaled down version of the Australian SKA for 
M = 120 antenna locations. Figures 3(a) and |3(b)] show the XY plane when 



populated with all the 120 antennas and after removal of 40 antennas and 80 
antennas respectively. 

Since the Z coordinates of the antennas also contribute to defining the UV 
plane, we provide the corresponding XZ plots along with the XY plots. Fig- 
3(c) and 3(d)| show the corresponding XZ planes. 



urcs 



Figures 3(e) and |3(f)| show the corresponding UV-planes. The hour angle 



used is 0° and the declination —30° 

Note that by our scheme, more antennas are removed towards the center 
of the log-spiral. Thus, it is more effective when N is large. Also observe 
that by making use of the same Minimum Variance Method, we can create a 
construction plan for all remainder phases of the telescope building which has 
a steadily improving UV-coverage. 

4.3 Problem 3: Mobile Antennas 

Statement: Suppose there are L pads (for mounting L antennas), but only 
M < L mobile antennas pre-placed on M pads. Suppose that an astronomer 
requests these M antennas for some observation. However, only N out of these 
M can be allocated (say, because of maintenance reasons, or some being needed 
for another experiment). Our problem is to choose N out of the M antenna 
pads that best approximate the desired UV pattern generated by fully functional 
pre-placed M antennas. Having done that, the antennas can be shifted to the 
recommended N pads. 

The parameters set for the previous problem in Subsection |4.2| carries through 
exactly. We did the experiment on Sub Millimeter Array data [5]. We assume 
the failure of 5 of the 24 antennas. Thus, in this setting, M = 24 and N = 19. 

shows the initial and final XY planes and |4(b)| shows the XZ plane, 
show the resultant UV distributions. Using the output of the algo- 



4(a) 



4(c) 



Figure 
Figure 

rithm, one can recommend which 19 pads can be used. 



5 Conclusions and Future Work 

In this paper we proposed and studied a Minimum Variance Method that can 
be used in many radio antenna placement scenarios. We applied the method 
in three different situations effectively (i) obtaining ab-initio smooth Gaussian 
UV distributions (ii) incremental construction of very large aperture synthesis 
arrays and (iii) achieving fault tolerance in mobile antennas. We also report 
experiments that indicate the usefulness of this method. In addition versatility, 
this method is quite efficient when compared to the brute force method or minor 
improvements thereof. 

In this paper, the only criterion that we considered for placing the antennas 
is image quality. It remains to come up with solutions that take into account 
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(a) XY plane-initial and 40 removed anten- (b) XY plane-initial and 80 removed anten- 
nas nas 
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(c) XZ plane-initial and 40 removed anten- (d) XZ plane-initial and 80 removed anten- 
nas nas 























/ 

t -J 










- -7 - 




588888 
--.HHHHH 










, :/ 
















F -;• 
'•</ •. 

• 4 















































































































-4 - 

-w * • 




588888 


























•V 


/ . 




































-588888 

































(e) The UV plane— initial and 40 removed an- (f ) The UV plane— initial and 80 removed an- 
tennas tcnnas 

Figure 3: Results for the proposed Australian SKA configuration, 
logistical factors like wire length minimization, roadway utilization and so on. 
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(a) The XY plane 



(b) The XZ plane 
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(c) The UV plane 



Figure 4: Results for SMA configuration (mobile antennas) 
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A Aperture Synthesis Primer 

It is a well-known fact in astronomy that the angular resolution of a telescope 
is proportional to X/D where A is the wavelength of the waves to be observed, 
and D is the diameter of the telescope aperture. In order to observe celestial 
bodies in the long wavelength ranges (1mm to 10m), we need antennas with 
diameters of the order of hundreds of kilometers. Obviously, such antennas 
are practically impossible to build, manoeuver and maintain. Fortunately, it 
has been shown that several antennas can be used instead of a single big one. 
Several radio telescopes have been constructed in the past based on this principle 
of aperture synthesis. These telescopes have helped astronomers in discovering 
various celestial bodies like quasars, pulsars etc. 

In this Appendix we quickly recap some concepts and terminology pertinent 
to this paper. The interested reader is referred to [3] and [5] for an excellent 
treatment of the subject. 

Let us first consider the case of two antennas. Consider a rectilinear coordi- 
nate system (u,v,w) such that the w-axis points in the direction that is to be 
the center of the synthesized field of view. Let So be the unit vector along the 
w-axis. 

Assigning a direction and treating the shortest distance between the two 
antennas as the magnitude, we can speak of a baseline vector [[] Let the baseline 
vector be b = uu + vv + ww. Let I, m and n be the direction cosines of an 
element in the celestial sphere. Let s = lu + mi + nw be a position vector of 
another element close by. Then, b.s = ul + vm + wn and b.So = w. 

Let A(l, m) be the response of the antennas corresponding to a baseline in 
the direction specified by (I, m). Further, let Aq be the response of the antennas 
at the beam center. Then, A(l, m) — is called the normalized antenna 

reception pattern. The spatial correlation function of the electric field at the 
antennas corresponding to the baseline b is defined as 



where 1(1, m) is the radio brightness in the (I, m) direction, assuming that both 
the antennas are identically polarized. 

Obtaining 1(1, m), our aim, is simplified to evaluating the inverse of a 2-D 
Fourier transform if w = : 



Therefore, we sample a plane perpendicular to w, the direction of our inter- 
est. Every pair of antennas generates a point on the UV-plane will be called a 
uv-point. 

3 In general, a long baseline improves the resolution of the array, while a short baseline 
implies a larger field of view. A larger number of different baselines implies higher sensitivity. 



V(u, v,w) = / / A(l,m)I(l,m)e 



2-Ki\ul+vm+w{s/ 1 — I 2 — m 2 — 1)] 



dldm 



J — OO J —OO 



yl — I 2 — m 2 
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Figure 5: The XY and UV planes 



In the case of a large number of antennas, a convenient coordinate system is 
used for the XY-plane [TT]. We choose the Z— direction to be along the north 
pole, the X— axis and the Y — axis lie on the equatorial plane with the X— axis 
pointing in the direction of the Greenwich meridian and the Y— axis pointing 
to the East, and the origin (0,0,0) lies at the center of the earth. Figure [5] shows 
the XY and UV planes in some detail. 

The following relates the latitudes and longitudes of a point on the earth to 
its X and Y coordinates in the system defined above. 

x = R e cos(e tat ) cos(e long ) 

Y = R e cos(6i a t)sm(0ion g ) 
Z = R e sm(0iat), 



where R e is the radius of the earth, and 6i a t and 9i ong are the latitude and 
longitude of the location respectively. 

Consider now two antennas placed at {X\,Y\,Z-\) and (X%, Y2, Z2). Let 
L Y = X 2 - X 1} L Z = Y 2 ~ Y 1 and L x = Z 2 - Z x . 

Then, the (u, v, w) coordinates are related to Lx, Ly , Lz as follows: 

sini? cosiJ \ / L x \ 

— sinScosH sin SsinH cos S Ly (1) 

cos ScosH — cos SsinH sin 5 I \ Lz I 

where H is the hour angle and S is the declination of the observation direction 
respectively, and A is the wavelength. While the curve traced out by a single 
baseline is coplanar, different baselines will lie on parallel planes, but at different 
elevations. The equation of a UV-curve is given by 



V 






-i( 






V w J 





I v — (L z /X) cos A 2 L\ + L\ 
sin 5 J A 2 
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The equation clearly shows that for given values of Lx , Ly and Lz , the path 
trace by a UV point is an ellipse. If the Lz = 0, the origin will be the center 
of the ellipse. For Lz ^ the ellipse center is shifted accordingly on v axis 
by (Lz/X) cos 5. For example, see Figure 6(a) and |6(b)] plotted for Australian 
proposed SKA antennas/stations. 




We finish by stating some useful definitions and remarks. The earth's axis 
determines half planes with the zenith and the source respectively. The angle 
between these two planes is called the hour angle. The angular distance of a 
source above or below the celestial equator is called its declination. 

For good image quality, it is desirable that we sample as many points on 
the UV-plane as possible. As the earth rotates, the UV-points trace out curves 
on the UV-plane, populating it. This technique of making use of the earth's 
rotation to fill up the sampling plane is called aperture synthesis. Naturally, we 
require that over time, the UV-plane is filled up as much as possible. 

It is a useful exercise to work out the curve traced by a UV-point as the 
earth rotates. In particular, one will observe that the curve traced for a source 
at declination is a straight line, that for 90° declination is a circle, while all 
intermediate sources trace out ellipses of increasing eccentricity. 

While different sources would require different UV distributions, it has been 
suggested that (i) in case of a dense UV plane, a Gaussian distribution is pre- 
ferred in the radial direction, and a uniform distribution is preferred across the 
azimuth; (ii) in case of a sparse UV plane, a uniform distribution over the plane 
is preferred [2] . 
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B Area Under the Gaussian 

_ x 2 

Consider a Gaussian y — ae ^ centered at x — 0, amplitude a and variance 
a 2 . Denote by A the area under the curve: A = a J_ e^^dx. Then, A 2 = 
a 2| r°° e '^ dx \ 2 . Or, A 2 - a 2 f°° f°° dady. 

I J — oo J ' J — oo J — oo a 

Writing it in polar coordinates, we get A 2 = a 2 f * d8 f^ e~2^ rdr. This 
evaluates to 2i:a 2 c 2 . Therefore, A = ^/2nac. 
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